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SOME EFFECTS OF FREQUENCY ON THE CONTRIBUTION OF A VERTICAL TAIL TO THE 
FREE AERODYNAMIC DAMPING OF A MODEL OSCILLATING IN YAW ' 

By John D. Bird, Lewis R. Fisher, and Sadie M. Hubbard 


SUMMARY 


SYMBOLS 


The directional damping and stability of a f uselage — vertical- 
tail model oscillating freely in yaw were measured at a Mach 
number of 0.1 4 and compared with the damping and stability 
obtained by consideration of the erffects of unsteady lift. The 
contribution of the vertical tail to the damping in yaw of the 
model agreed within the limits of experimental accuracy with 
predictions of the approxinoate finite-aspect-ratio unsteady-lift 
theories for the range of freguencies tested, and the two- 
dimensional unsteady-lift theory predicted a much greater loss 
in damping with reduction in freguency than was shown by 
experiment or by the approximate finite-aspect-ratio unsteady-lift 
theories. For freguencies comparable to those of lateral airplane 
motions, both the unsteady-lift theories and the experimental 
results indicated that the contribution of the vertical tail to the 
directional stability of the model was relatively independent of 
freguency, 

INTRODUCTION 

The advent of high-speed airplanes of high relative density 
has focused attention on certain problems associated with 
the dynamic stability of airplanes which, because of previous 
unimportance, have heretofore been neglected. Among these 
problems is the effect of the periodicity of the airplane motion 
on the effective values of the various stability derivatives. 
Reference 1 and, more recently, references 2 and 3 have 
indicated the possibility of sizeable effects from this problem. 

An appreciable amount of theoretical work on these 
unsteady-lift effects exists at present; however, only a small 
amount of experimental substantiation of the results is 
available for the low-frequency range of oscillation. As a 
result, a program has been undertaken in the Langley sta- 
bility tunnel to determine the effects of such variables as 
frequency and amplitude of motion on the contribution of 
the various airplane components to the stability derivatives 
of present-day airplane configurations. The work reported 
herein covers that phase of the investigation which considers 
frequency effects on the directional damping and stability 
of a model undergoing a freely damped oscillatory yawing 
motion. The effects of vertical-tail aspect ratio and com- 
pressibility as predicted by the theoretical treatments are 
discussed in relation to the experimental stability charac- 
teristics obtained by the free-oscillation and by the curved- 
fiow procedures. 


The data are referred to the system of stability axes and 
are presented in the form of standard NACA coefficients of 
forces and moments about the quarter-chord point of the 
mean aerodynamic chord of the normal wing location of the 
model tested. (See fig. 1.) The coefficients and symbols 
used herein are defined as follows: 

Cl lift coefficient, LjgSw 

Cd drag coefficient, DjgSy, 

Cy lateral-force coefficient, YlgSw 

Cl rolling-moment coefficient, L' jgSJiy^ 

Cm pitching-moment coefficient, MjgS^jfCy, 

Cn yawing-moment coefficient, NlgSJ)w 




Section A- A 

Figure 1. — System of stability axes. Arrows indicate positive forces, 
moments, and angular displacements. Yaw reference generally 
chosen to coincide with initial relative wind. 


1 Supersedes NACA TN 2657, “Some Effects of Frequency on the Contribution of a Vertical Tail to the Free Aerodynamic Damping of a Model Oscillating in Yaw’* by John D. Bird, 
Lewis R. Fisher, and Sadie M. Hubbard, 1952, 
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C{k)=F+ia 

F,G 


F=F+iO) 

Q=H+ij] 


Z,M 


lift, lb 
drag, Ib 
lateral force, lb 
rolling moment, ft~lb 
pitching moment, ft-lb 
yawing moment, ft-lb 

dynamic pressure, p F^, Ib/sq ft 


mass density of air, slugs/cu ft 
free-stream velocity, ft/sec 
area, sq ft 
span, ft 
chord, ft 

mean tail chord, ft 


2 

mean aerodynamic chord, ^ J cHy, 


ft 


spanwise distance from the plane 
symmetry, ft 
aspect ratio, ¥jS 
angle of attack, deg 
angle of sideslip, radians 


of 


time, sec 


3 ^awing-velocity parameter referred to 
semispan of wing 
angle of yaw, radians 
yawing velocity, dypldtj radians/sec 
yawing acceleration, d'^ypldt^j radians/sec^ 
angular velocity, radians/sec 
frequency, cps 

reduced-frequency parameter referred to 
semichord of vertical tail, cocj2F 
mechanical spring constant, ft-lb/radian 
3 ' awing moment of inertia, ft-lb-sec^ 
time to damp to one-half amplitude, sec 
daniping constant 
period of oscillation, sec 

nondimensional tail length referred to semi- 


chord of vertical tail, ^ 

Ctl2 

distance from origin of axes to midchord 
point of vertical tail, ft 


circulation functions used 1)3" Theodorsen 
(ref. 6) 

finite-span functions used by Biot and 
Boehnlein (ref. 9) 

finite-span correction used by Reissner and 
Stevens (ref. 8) 

compressible-flow functions used in place of 
and G (refs. 10 and 11) 


A,B functions of k defined in text 

^ rb 


1>JV‘ 




bCv 


rb 

W 


r _ ^Cn 

dp 


On, 


d 

2F 


Cy = 


bCy 


A^i= 


N2= 


dp 

di^ 

d^ 


AT 

AT 

N-p = —r 


Subscripts: 
w wing 

t vertical tail 

0 amplitude 

f due to friction 

APPARATUS 

MODEL AND OSCILLATING STRUT 

The model used in this investigation and shown in figure 
2 had separable wing, fuselage, and tail surfaces and was 
constructed almost entirel 3 ' of balsa wood in order to mini- 
mize the mass and make tlie natural frequency of the model 
on its mounting as high as possible. The configuration is 
representative of that of several present-day high-speed 
airplanes. The aspect ratio At and mean chord Ct of the 
vertical tail were 2 and 4.9 inches, respectively. 

The model was mounted at the quarter chord of the wing 
mean aerodynamic chord on a 1 -inch-diameter rod which 
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Torque rod " 


A Recorder 
B Lens 
C Film 
D Timing light 



Figure 3. — Diagram of test equipmeut. 



Figure 2. — Drawing of model tested. All dimensions are in inches. 

rotated with respect to a circular hollow strut and which 
was connected to it by means of flexure pivots made of 
Swedish steel. Flexure pivots were used rather than bear- 
ings in order to minimize the friction in the system. The 
pivots were alined with the center of the rod on which the 
model was mounted. The entire assembly, which is shown 
in figures 3 and 4, was fastened to the ceiling of the tunnel 
test section and covered by a streamlined fairing. 

The inner rod or oscillating strut, which was permitted to 
rotate and upon which the model was mounted, extended 
through the ceiling of the tunnel test section and was capped 
by a horizontal crosspiece. The period of oscillation of the 
model was controlled by varying the moment of inertia 


about the vertical axis of the system. This variation was 
accomplished by using crosspieces of different lengths and 
by clamping weights at the ends of the crosspieces. 

RECORDING DEVICE 

A continuous record of the displacement in yaw of the 
model was made on photographic film by means of an optical 
recording system. A mirror clamped to the section of the 
oscillating strut which extended outside the tuimcl reflected 
a beam of light from a point light source within the recorder. 
The beam was focused on a roll of photographic film liy 
means of an 8-inch-focal-length double-concave lens. The 
film roll was driven at constant speed by gearing it to an 
electric motor. A 10-cyclc-pcr-second neon timing light 
incorporated in the recorder exposed timing lines on the film 
in order that time as well as model displacement could be 
measured accurately on the film record. 

TESTS 

OSCILLATION TESTS 

The free-damping tests were conducted with the model in 
an inverted position, in order to minimize strut interference 
effects, and at zero angle of attack with respect to the 
fuselage reference line. The tests were conducted on the 
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Figure 4. — Test model mounted on free-oscillation strut in Langley 
stability tunnel. 


combination of fuselage and tail surface and on the fuselage 
alone because preliminary tests had shown that the presence 
of the wing had no apparent effect on the damping at zero 
angle of attack. The range of moments of inertia and fre- 
quencies covered are shown together with the test results 
in table I. Results of runs made with wind on and off are 
given in order to illustrate the magnitude of the tares 
applied. 


^Yllen tests were made of the fuselage alone, a torque rod 
was em()lo 3 ^ed to provide a stable restoring moment to the 
fuselage which was equal to the aerodynamic moment pro- 
vided b}^ the vertical tail when attached. The torque rod 
was clamped between the extended oscillating strut and a 
building structural member as shown in figure 3. The 
sjiring constant of the S 3 "stem was G.8 foot-pounds per radian 
for the flexure pivots alone and 32.8 foot-pounds per radian 
for the coml)ination of the flexure pivots and the torque rod. 

FORCE TESTS 

The stead 3 ^-state stability charactei-istics of the model 
were determined by means of wind-tunnel force tests through 
an angle-of-attack range from 0° to 8° for the fuselage alone, 
the fuselage and tail surfaces, and the complete configura- 
tion. The model was mounted on a single-strut su])port and 
all forces and moments were measured by a six-component 
balance system. The steady-state sideslipping derivatives 
were determined through a sideslip range of ±5°, and the 
yawing <lerivatives were detennined at values of the 3 Tiwing- 
velocity parameter Thj2V of 0, — 0.0291, — 0.0016, and 
— 0.0810 by standard Langle 3 ' stability-tunnel curved-llow 
testing procedures. 

All tests were conducted in the 6- by^ 6-foot test section of 
the Langley stability tunnel at a dymamic pressure of 24.9 
pounds per square foot, a Alach number of about 0.14, and 
a Reymolds number of 442,000 based on the wing mean 
aerodymamic chord. Standard jet-boundaiy corrections 
have been applied to the angle-of-attack, the drag-coeflicient 
and the rolling-moment-coefficient data obtained from the 
force tests. In addition, the steady^-state y^awing data have 
been corrected for the effects of the lateral static-pressure 
gradient peculiar to the curved-flow testing j)rocedure. 
Jet-boimdaiy corrections were not applied to the oscillation 
tests because of their insignificance for the test conditions 
involved. See reference 4. 

REDUCTION OF OSCILLATION TEST DATA 

From the continuous film record taken of the motion of the 
model after an initial displacement, the amplitudes of the 


TABLE I 


EXPERIMENTAL RESULTS 


Configuration 

p, 

sec 

k 

h/2, 

sec 

h, 

ft-lb-sec 2 


(Cn.-C„j) 



Wind 

off 

Wind 

on 

Wind 

off 

Wind 

on 

Total 

Friction 

Aero- 

dynamic 


13. 21 

7. 15 

0. 0012 

190. 08 

85. 80 

32. 89 

0. 610 

0. 275 

0. 335 

0. 204 


4. 08 

2. 19 

. 0040 

49. 70 

11. 47 

3. 09 

. 429 

. 099 

. 330 

. 204 

Fuselage and tail 

3. 96 

2. 25 

. 0039 

54. 65 

10. 35 

3. 09 

. 475 

. 090 

. 385 

. 190 


3. 96 

2. 19 

. 0040 

54. 65 

11. 39 

3. 09 

. 432 

. 090 

. 342 

. 204 


. 75 

. 77 

. 0114 

34. 36 

1. 72 

. 38 

. 351 

. 018 

.333 

. 213 


. 22 

. 30 

. 0293 

3. 40 

. 28 

. 06 

. 340 

. 028 

. 312 

. 203 


. 29 

. 27 

. 0321 

13. 80 

3. 21 

. 06 

. 030 

. 007 

. 023 

-. 058 

Fuselage alone 

. 75 

. 74 

. 0119 

34. 36 

13. 77 

. 38 

. 044 

. 018 

. 026 

015 
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successive cycles were measured and plotted to a logarithmic 
scale against time. Inasmuch as the damping is of a loga- 
rithmic nature, the resulting plot is a straightline from which 
may be read tlie time for the motion to damp to one-half 
amplitude (/ 1 / 2 ). 

For a system having a single degi’ee of freedom in yaw 
suclv that for the stability-axis system the angle of sideslip /3 
is the negative of the angle of 3 ^aw xj/j the equation of motion 
can be expressed as 


where and N^p refer to the mechanical friction of the 
system and the flexure-pivot spring constant, respectively. 
The solution of equation (I) can be written as 


where 


yp=ze~^^{A sin 27rft+B cos 27 t/0 


/ n,-n, I 

rVi-iVi+ViY 

' Iz-N-i \ 

L 2{h-N-i) J 


m = 




( 2 ) 


The damping constant may be expressed in terms of the time 
to damp to one-half amplitude as 


0.693 
m=— 

H/2 

The directional stability and damping of the system may 
then be expressed in nondimensional form in a manner similar 
to that of reference 5 as 




c c ‘ 1-^i 

gSwbt/ \ L^ 1/2 {t'l/^windoffA ^1/2 j 


( 4 ) 


The terms Cn- and Cn^ are necessarily in combination with 
Cn^ and respectively, because of the nature of the motion 
being considered. For the frequencies employed in these 

tests, the factor should be small in comparison 


with Orig- 
in determining (ti, 2 )wind ofh which corrects for the damping 
due to friction in the system, the model was oscillated at 
zero wind velocity with an equivalent mass in place of the 
vertical tail and with a torque rod of sufficient strength to 
produce the proper frequency inserted in the system. The 

influence of the term is very small compared with that of 

Q/2 

the inertia term in equation (4) and was neglected in these 
tests. 

The damping derivatives were obtained from data similar 
to the sample records of model motion and plots of logarithm 
of amplitude against time given as figure 5. A small residual 
yawing motion of the model attributable to turbulence in the 
tunnel airstream necessitated determining the time to damp 


to one-half amplitude by the slope of the amplitude envelo])e 
at the. largest angles of yaw in order to minimize the error in 
determining the damping derivatives. This source is 
believed to have resulted in no more than a 5-percent 
uncertainty in the results. 

ANALYTICAL CONSIDERATIONS 

Various unsteady-lift theories were employed to calculate 
the contributions of the vertical tail to the directional 
stability and damping of the model for comparison with the 
experimental results. These theories included the incom- 
pressible cases covered by Theodorsen (ref. 6) for two 
dimensions and by Jones (ref. 7), Reissner and Stevens (ref. 
8), and Biot and Boclmlein (ref. 9) with various degrees of 
approximation for finite span. The two-dimensional work 
including the subsonic effects of compressibility as given by 
Possio and Frazer and Skan is also employed. A short 
smnmary indicating procedures for obtaining the necessary de- 
rivatives from these theories is included here for completeness. 

TWO-DIMENSIONAL DERIVATIVES 

The work of Theodorsen on the derivation of the expres- 
sions for the unsteady lift and moment of a two-dimensional 
surface undergoing sinusoidal oscillations in an incomjires- 
sible fluid is given in reference 6. From this source, the 
moment per unit span about the airplane center of gravity 
of a vertical tail mounted on an airplane which is performing 
sideslipping oscillations can be shown to be 

t7(t)/3] (5) 

A similar expression for an airplane performing }^awing 
oscillations is 

w)— I + j (!+«■) i,] (6) 

An ajipendLx to this report gives a discussion of the con- 
version of these results from reference 6. 



(Q) 



(a) Film record of yawing oscillation. 

(bj Logarithmic plots of amplitude variation with time. 
Figure 5. — Representative film record and resulting logarithmic plots. 
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Tlie Theodorscn function C{k) is a complex circulation 
function of the reduced frequency and in terms of its real 
and imaginary parts is written 


C{k)=F+iG 

The functions F and G may be computed from expressions 
given in reference G in terms of Bessel functions of the first 
and second Idnds of orders zero and one and of argument 
k. In stability practice equations (5) and (6) are usually 
interpreted in terms of derivatives as follows: 

N{0) = + C„. ^ 

mi) = H Cn, 'i j 

If the sinusoidal motions 



are emplo 3 "ed for a substitution into equations (5), (6), and 
(7), the real and imaginaiy terms may be equated and the 
derivatives Cn„ and Cnj. determined. The results of 
such an analysis would be as follows: 




St Cl 


( k ) ' [I +(“+0 f ] 


(8) 


The velocity terms contain only the F term of the Theodorscn 
function, whereas the acceleration terms contain the G term 
divided b}^ k. The corresponding side-force flerivatives are 
derived in the appendix. For the free-oscillation tests 
conducted for this report, the motion is such that the angle 
of sideslip is the negative of the angle of yaw. The damping 
derivative measured therefore is the difference between the 
derivatives C'nr determined from equation (8). 

Tliis fact Avas illustrated previousl}^ with respect to the 
reduction of the oscillation test data. In a like manner, the 
directional stability parameter measured can be shown to be 
a combination of the and Cn; terms, specifically 


Cn^+Cn,k^ 



Combination of the results of equations (8) to obtain the 
factors calculated for comparison with the oscillation tests 
gives 


r —r — — 

2\bJ s^k 


( 9 ) 


wlicrc' 

^=-(a-|) («+|) 2kF 

2k6 

FINITE-SPAN DERIVATIVES 

The aero(hmamic-span effect is considered in reference 7 
for wings of as])ect ratios of G and 3 by correcting the aero- 
dvnamic inertia and the angle of attack of the infinite-span 
surface. An approximation employed in this anatysis 
concerning motions of long wave length may make the 
results subject to question for values of k near zero. The 
reference presents expressions for calculating the finite-span 
values of the Theodorscn function C(k). The use of these 
F and G circulation functions in equations (8) and (9) 
results in values of the various derivatives considered for 
as])ect ratios of G and 3. 

The aerod^mamic-span effect is considered for wings of 
arbitrary aspect ratio in reference 8. The three-dimensional 
effect of the finite span ma}^ be obtained b^^ adding a correc- 
tion term a to the basic Theodorscn two-dimensional func- 
tion. The finite-span function C(k)~{-o- is emplo}"ed in a 
rather lengtly and tedious computation in order to deter- 
mine the components of the unsteady lift for each of five 
points along the span of the oscillating surface. The various 
stability derivatives may then be determined b}- a graphical 
integration of the in-phase and out-of-phase components. 
An approximation made in this analysis is of such a nature 
as to make the method primarily suited to the higher aspect 
ratios. A practical lower limit of applicability is not known. 

A one-point approximation to the span effect for the cal- 
culation of the unsteady lift may be employed rather than 
the five-point procedure. Comparative calculations made 
for the present investigation have shown the one-point 
approximation to ^deld results for low aspect ratios which 
are essentiall}’^ the same as those calculated b^- the longer 
procedure. 

Biot and Boelmlein (ref. 9) considered the aerodAuiamic- 
spaii effect on unsteady lift by an approach that yields a 
closed form for a one-point approximation to the span effect. 
The finite-span circulation functions used aj*e 

Q=H+iJ 

For infinite aspect ratio 

p=Q=.C{k) 

where C{k) is the two-dimensional Theodorsen function. 
The functions P and may be calculated for the argument 
k by methods given in reference 9, and equations (8). and 
(9) are again used to calculate the various stability deriva- 
tives. The factors A and B which appear in equations (9) 
must, however, be written in terms of P and Q as follows: 

A=(2a-l)akG+(^a~^^ kj+(^a^+^^ k^+2aF+H 
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B={2a—\)akF-\-(^a—^ kll—(^a—^ k—2aG—J 

EFFECT OF COMPRESSIBILITY 

Subsonic com])ressibility effects on tlie unsteady dam])ing 
in yaw and directional stability of a two-dimensional oscil- 
lating lifting surface can be considered tlirough the use of 
reference 10 which is an extension of the calculations of 
Possio (ref. 11). Tlie infinite-span compressible-flow values 
for the stability derivatives can be obtained by the use of 
equations (8) and (9) for Mach nuinl)ers 0, 0.5, and 0.7. In 
the notation of Prazer and Skan (ref. 10), however, 

.1 = 4 0 ) ( + z,-| Z.)+ A/. 

/y=4A/,-l(l+a)(A/2+Z4-|Z,)+A/2 

The correspondence between the notation of Tlieodorsen 
and that of Frazer and Skan is shown in reference 12. The 
values of the Z and M functions are tabulated by Frazer and 
Skan for Mach numbers 0, 0.5, and 0.7 and a significant 
range of reduced frequencies. Reissner has recently ex- 
tended his work on oscillating wings of finite span to include, 
to the same degree of approximation as for the incompressible 
case, the subsonic effects of coin])ressibility. (See ref. 13.) 

RESULTS AND DISCUSSION 

PRESENTATION OF RESULTS 

The experimimtal values of the longitudinal, directional, 
and yawing characteristics of the complete model, the model 
with wing removed, and the fuselage alone as determined by 
force tests are given in figure 6. The experimental values of 
the damping in yaw of the fuselage and tail assembl}" 
(designated FT) and the fuselage alone (designated F) 
obtained by the freely damped oscillation technique for 
several frequencies of oscillation are shown in figures 7 and 8. 
The frequency and aspect-ratio effects on the damping in 
yaw contributed by the vertical tail as calculated from the 
imsteady-lift theories discussed previously are also shown in 
these figures. Figure 9 gives a comparison of the relative 
importance of Cnj. and Cn^ on the damping in yaw of the 
model. Figures 10 to 12 are similar figures summarizing 
frequency and aspect-ratio effects on directional stability. 
Figures 13 and 14 give the computed effects of compressi- 
bility on the vertical-tail contribution to the damping in yaw 
and directional stability of the test model. 

DISCUSSION 

The variation of the experimentally determined damping 
in yaw of the test model with frequency is in reasonably 
good agreement for the range covered by the tests with the 
various approximate theoretical treatments which consider 
the effect of finite span. (See figs. 7 and 8.) The contribu- 
tion of the fuselage to the damping in A^aw is small for the 
range of test frequencies. A two-dimensional approxima- 
tion to the damping as is frequently employed for flutter 
work appears to be remarkably good for and for aspect 


ratios of 6 and greater (fig. 8). A large effect of aspect ratio, 
however, is shown for aspect ratios below 6 for the entire 
frequency range and for all aspect ratios for fi*equencies less 
than ^=0.1. The two-dimensional result indicates a loss in 
dam])ing at low frequencies much greater than the finite- 
span result and, in fact, shows a reversal in sign of the 
damping at sufficiently low frequencies. A number of papers 
dealing with unstead}^ lift have noted this chai’acteristic of 
the two-dimensional aerod 3 mamic forces. (See refs. 1 to 3 
and 14.) 

An examination of the values of the derivatives Cn^ and 
as computed from reference 6 for infinite aspect ratio 
and from reference 7 for aspect ratios of 3 and 6 (fig. 9) 
indicates that the derivative is criticall.v dependent on 
reduced frequenc.v for the low range of reduced frequencies 
and infinite aspect ratio. The derivative Cn^ is only slightly 
dependent on the reduced frequenc}^ for infinite aspect ratio 
and even less so for the lower aspect ratios. The dependence 
of on reduced frequency is responsible for the decrease 
in damping shown for the infinite-aspect-ratio case as the 
frequency is reduced (figs. 7 and 8). 

The aspect-ratio-2 results shown in figures 7 and 10 were 
obtained from the theoiy of Jones (ref. 7) for aspect ratio 
3 b\" use of a correction factor which was the ratio of the 
lifting-line-theorv lift-curve slopes for aspect ratios of 2 and 
3. This correction is only approximate and was made in 
order to place the results on a comparable basis. The varia- 
tion of Cn^ — Cn-^ with reduced frequenc 3 ", shown in figure 7, 
indicates no appreciable effect of frequency to zero frequency'. 
In fact, the result of Jones for aspect ratio 3 reduces to a 
finite value at k=0 given by^ 

y-C„. = -|(^^y|^(1.2a^-0.370a + 0.515) 

which is not the case for the two-dimensional result that 
approaches infinity as k approaches zero. This result is not 
expected to be veiy accurate close to k=Q, however, because 
of the nature of the assumptions made in its derivation as 
mentioned previoush^ 

When the finite-span effects are considered, it appears that 
these effects are especially^ significant with regard to the 
damping in yaw of this model for aspect ratios smaller than 
6. (See fig. 8.) The theory also indicates a significant 
decrease in the damping at low frequencies when the aspect 
ratio was increased from 20 to » . 

The experimental values of damping in yaw obtained by 
free oscillation agree closely with the values of Cn^ obtained 
by- the standard curved-flow testing procedure used in the 
Langley stability tunnel. The curved-flow technique results 
in the ^=0 value of Cn^. No unsteady effects, however, 
are obtained by^ this technique. 

The unsteady-lift theories indicate that, at frequencies of 
oscillation from A:=0 to A'=0.1, the frequency' has little 
effect on the directional stability contributed by the vertical 
tail but that high frequencies of oscillation can cause the 
directional stability contributed by the vertical tail to de- 
crease greatly and eventually to change sign (figs. 10 and 11). 


206367 — 54 ^ 
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Fioure 6. — The steady-state stability characteristics of the moflel tested. 
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Figure 7. — The effect of frequency on the damping in yaw for a low range of the reduced frequency. Comparison of several methods for cal- 
culating frequency effect on the damping in yaw. 
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Fir.uRE 8. — The effect, of aspect ratio on the damping in yaw as calculated by the method of Biot and Boehnlein for a range of reduced frequencies. 
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(b) Range of reduced frequencies. 

Figure 9. — The effect of reduced frequency on the stability derivatives 

and 

Changes in aspect ratio, for low aspect ratios, have a large 
effect on the directional stability; for high aspect ratios, 
changes in aspect ratio are less significant. The experimental 
values of directional stability obtained from the free-damping 
tests remained constant for the range of frequencies investi- 
gated, although the magnitudes are somewhat smaller than 
those predicted for the tail contribution by the calculation 
methods for aspect ratio 2. This result is not unexpected 
because all the finite-span unsteady-lift theories employed 
maintain some resemblance to the Prandtl lifting-line 
theory which, of course, yields values of the lift-curve slope 
for low aspect ratios which are unrealistically high. The 
values of directional stability obtained by both the free- 
oscillation and the static-test procedures were about the 
same magnitude. 

Figure 12 gives a comparison of the contribution of a 
vertical tail to the stability derivatives Cn^ and (7^; as 


computed from references 6 and 7. These derivatives were 
combined in the proportion 

for comparison with the results of the oscillation tests. 
Although Cn'^ is of appreciable magnitude itself, the factor 

small relative to for the range of 

frequencies around fc=0.01 and thus should not greatly 
affect the directional stabilit}^ of the usual airplane. 

The computed two-dimensional Mach number effects on 
the damping in yaw and directional stability of the tail 
surface of the test model are shown in figures 13 and 14. 
These results were included to illustrate the powerful effect 
of compressibility on the damping contributed by the tail 
surface. A complete coverage of one-degree-of-freedorn 
flutter by use of these aerodynamic forces which may be 
obtained from references 10 to 12 is given in reference 14. 
Compressibility effects decrease the tail-surface damping at 
low frequencies and increase the damping at higher frequen- 
cies. The directional stability increases at all frequencies 
with increase in Mach number with the largest increases 
occurring in the higb-frequency range of oscillations. Refer- 
ence 15 which presents the two-dimensional oscillating 
aerodynamic forces at sonic speed indicates larger deviations 
from the results for incompressible flow than for the subsonic 
results given here. 

CONCLUSIONS 

The following conclusions were drawn from the results of 
an investigation made at a Mach number of 0.14 of the 
effects of frequency and aspect ratio on the directional 
damping and stability of a fuselage and vertical tail oscil- 
lating freely in yaw: 

1. The contribution of the vertical tail to the damping in 
yaw of the model agi^eed within the limits of experimental 
accuracy with predictions of the approximate finite-aspect- 
ratio unsteady-lift theories for the range of frequencies 
tested. 

2. Two-dimensional unsteady-lift theory predicted a much 
greater loss in damping with reduction in frequency than was 
shown by experiment or by the approximate finite-aspect- 
ratio unsteady-lift theories. 

3. Both the unsteady-lift theories and the experimental 
results indicated that the contribution of the vertical tail to 
the directional stability of the model was relatively independ- 
ent of frequency for frequencies comparable to those of 
lateral airplane motions. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 26^ 1951, 
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Fioube 10.— The effect of frequency on the directional stability for a low range of the reduced frequency. Comparison of several methods for 

calculating frequency effect on the directional stability. 
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Figure 11. — The effect of aspect ratio on the directional stability as calculated by the method of Biot and Boehnlein 

for a range of reduced frequencies. 
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(a) Lower range of reduced frequencies. 

(b) Range of reduced frequencies. 

Figure 12. — The effect of reduced frequency on the stability deriva- 
tives Cnp and Cn-> 



(a) Lower range of reduced frequencies. 

(b) Range of reduced frequencies. 

Figure 13. — Tlio effect of compressibility on the damping in yaw. 
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(a) Lower range of reduced frequencies. 

(b) Range of reduced frequencies. 

Figure 14. — The effect of compressibility on the directional stability. 


APPENDIX 

CONVERSION OF UNSTEADY FORCES AND MOMENTS FROM FLUTTER TO STABILITY AXES 


The conversion of the contribution of the vertical tail to 
the yawing moment of an airplane in terms of unsteady-lift 
functions and of the Theodorsen flutter axes to the stability 
axes which are fixed in the body is included here as a matter 
of general interest. Reference 16 gives considerable atten- 
tion to such axis conversions. The following sketches indi- 
cate the paths and attitudes of airplanes performing sinusoidal 
sideslipping and yawing oscillations within the limits of the 
stability-axis system: 




Yawing, 

The Theodorsen coordinates h and a and the stability-axis 
coordinates /3 and \f/ are shown with respect to the paths 
followed by the airplanes in sideslipping and yawing flight. 
Ill sideslipping flight 


Sideslipping, j3=jSoe‘"‘ 


a=0 


16 

and in yawing flight 
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— -^h = a 


From reference 6 the contribution of the vertical tail to the 
yawing moment of an airplane in terms of the flutter coordi- 
nates is 

,V , (o+j) cm [F«+i+| (j-®) i] 

If this moment is for a sideslipping motion or a yawing 
motion in the sense given by stability axes fixed in the body 
and if the following terms are substituted for the sideslipping 
case 

a= 0 
a= 0 
a = 0 

h=-Vfi 

h=-V^ 

and for the yawing case 

h=-^V^ 

h=-V^ 

a=\J/ 
a=\l/ 
a = \l/ 

the following results are obtained for the sideslipping case 

iV(^)=-^ 4)+^^ (a+i) C{k){-v^) 


awCiV 


and for the yawing case 

(i+«=) i- 

^ ft 

which reduce for the sideslipping case to 

N{0) = -irqc,^ Cm~\ 

and for the yawing case to 

m*)—l ,C.’ {[i-ew(i-.=)]i+^(i+a>)|, 




^ + 


sideslipping and yawing flight. Reference 6 gives the side 
force of the vertical tail in terms of the flutter axes as 

Y{a,h)=^ (yxa+x^-’^)+ 

TpVc,C{k) [ya+A+§ G““) “] 

wliich, when the following terms are substituted for the 
sideslipping case 

h=-Vfi 
h = -V^ 


and for the yawing case 


a=0 
a = Q 
a = 0 

h=-VyP 

h=-V^ 

a=}p 
a = \l/ 

OL = \p 


becomes for pure sideslipping and yawing flight 

r W= 3C.^ I \_2C{k) (i-a) I a |,] 


and 


Substituting the motions, in the manner done in the text, 
for the yawing moment, 

into these expressions and into the expressions in terms of 
stability derivatives which are 

Y(0) = Cr^^qS^+Cy.fi 

gives the following relationships for the tail contribution to 
the side-force derivatives when the real and imaginary parts 
of these equations are equated: 

(7. 


s„- 


Cy= — 7T 


^ C 

S* h 




— it (it) T G““)] 


Similar expressions may be obtained for the side force in pure 
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